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Abstract. The worldline formalism has been widely used to compute physi- 
cal quantities in quantum field theory. However, applications of this formalism 
to quantum fields in the presence of boundaries have been studied only re- 
cently. In this article we show how to compute in the worldline approach the 
heat kernel expansion for a scalar field with boundary conditions of Robin 
type. In order to describe how this mechanism works, we compute the contri- 
butions due to the boundary conditions to the coefficients Ai, A 3 / 2 and ^2 
of the heat kernel expansion of a scalar field on the positive real line. 



1. Introduction 

Many problems in quantum field theory can be studied in the worldline for- 
malism. This approach provides an intuitive and, in many cases, more efficient 
method to carry out computations in the I-loop approximation (see [21 13 E] and 
references therein.) 

This formalism makes use of the connection between some quantities in quan- 
tum field theory like the effective action, and the operator e~ TA corresponding 
to some relevant differential operator A. The first important step in this approach 
is to identify a classical particle whose Hamiltonian, after first quantization, is 
this differential operator A. If this identification is found, then the heat kernel 
corresponds to the evolution operator of the particle in Euclidean time T. The 
kernel (y\e~ TA \x) of this operator e~ TA is a quantum mechanical transition am- 
plitude and therefore can be also computed in the path integral formulation of 
first quantization. 

It is important to mention that many physical quantities of the quantum 
field theory are determined by the small-T asymptotic expansion of the hcat- 
kernel (y\e~ TA \x) . In this sense, the path integral computation of this transition 
amplitude results particularly appropriate because the small-T expansion comes 
out in a natural way. By this procedure the worldline formalism has been proved 
useful in the study of quantum fields of different spins and their coupling to external 
gravity. 

The purpose of this article is to describe how this formalism can also be 
applied to quantum field theory on manifolds with boundaries. In order to do this, 
one has to deal with the path integral formulation of the quantum mechanics of a 
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point particle on a manifold with boundaries. The main difficulty is to implement 
the boundary conditions imposed on the quantum fields, i.e. on the domain of the 
differential operator A, in the path integral calculation of the quantum mechanical 
transition amplitude (y\e~ TA \x) . This problem has been solved in [I] for Dirichlet 
and Neumann boundary conditions by the method of images. In this article we 
show how to implement Robin boundary conditions in the path integral calculation 
and we compute the contributions of the boundary conditions to the coefficients 
, A 3 / 2 and A2 of the small-T asymptotic expansion of the heat kernel for a 
scalar field on the positive half-line. For this setting, Robin boundary conditions 
are the most general boundary conditions determined by the selfadjointness of the 
differential operator A. 

The method of images is not straightforwardly applicable to Robin boundary 
conditions. However, as it is shown in the appendix (see also [S]), one can imple- 
ment Robin boundary conditions on the heat kernel by using the image construc- 
tion for Neumann boundary conditions supplemented by a suitable delta-function 
potential at the boundary. This delta-function then constrains the paths which 
contribute to the path integral computation of the transition amplitude. 

In the next section, we describe the worldline method by considering the 
heat-kernel trace for a scalar field on a manifold without boundaries. In section [3] 
we apply the worldline approach to a scalar field on the positive half-line under 
Robin boundary conditions at the origin and we compute the contributions of the 
boundary to the first coefficients of the heat kernel expansion. In section[4]we draw 
our conclusions and comment on possible generalizations of our technique. 



In order to describe the worldline approach to quantum field theories let us 
consider the 1-loop effective action of a scalar field <j>(x) with a self-interaction 
U((f>(x)) on a flat manifold A4. The classical action I[(f>(x)} of the scalar field is 
given by 



2. Heat-kernel expansion without boundaries 



(2.1) 




and the 1-loop effective action can be represented as 



T[(j)(x)} = I[<t>{x)] -hlogDet- 1 / 2 {-A + U"((f>(x))} 



(2.2) 




The integrand in the last expression is the heat kernel 



(2.3) (y\e- T i- A+u "^}\x) 
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of the Schrodinger operator 

(2.4) - A + U"(<p(x)} 

evaluated at the diagonal x = y. Inspection of expression (|2.2p shows that the 
small-T asymptotic expansion of this heat-kernel contains information about the 
ultraviolet divergences of the theory. 

In the worldline approach, the kernel (|2.3[) is regarded as the quantum me- 
chanical transition amplitude in Euclidean time T between x and y of a point 
particle whose Hamiltonian is given by (12. 4|) . In fact, the Schrodinger operator 
(|2.4j) can be regarded as the Hamiltonian, after first quantization, of a point par- 
ticle whose Euclidean classical action is 

(2.5) S[x(t)]=£dt [±x 2 (t) + V(x(t)) 
where 

(2.6) V(x) := U"{(j){x)) 

Therefore, we can represent the transition amplitude (12. 3|) in the path integral 
approach 

(2.7) (y\e- T ^ A+v ^\x) = / Vxit) e~ & dt U* 2 +n*«)} 

With this procedure we can compute some quantities in quantum field theory in 
terms of quantum- mechanical path integrals. An essential step in this description is 
the identification of the differential operator appearing in the first line of expression 
(|2.2p with the quantized Hamiltonian of a point particle. 

Next, we show how the path integral in expression (|2.7|) turns out to be con- 
venient for the computation of the small-T asymptotic expansion of the transition 
amplitude. In order to do that, we rescale the Euclidean time variable t by in- 
troducing r := t/T and we replace the integration over trajectories x(t) by an 
integration over the "quantum fluctuations" y(r) := x{Tt) — (y — x)t — x 

(y\e- T t- A+v ^\x) 



(2.8) 



e 4T 



'v(o)=o 

X e - T /o dTV{{y-x)T+x+y{T)) 

I -T^drV((y-x)T+x+y(T)) 



e it 



E 

n=0 



(-T)" f 1 



dr„ x 



l[V({y-x)T k +x + y(T k )) 



\k=l 
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Notice that we end up with a perturbative calculation in (0+l)-ficld theory where 
the number of loops is related to the power of T and the propagator is proportional 
to T. 

For our purposes it suffices to consider the trace of the heat kernel 

Tre -T{-A+V(x)} _ J dx (z-T ft dTV(x+y(r)) S 

(2.9) = Y, -^- I dx I dT i-' - / dT » ( \ \V(x + y( Tk )) 



n=0 n - JM Jo -„ u=1 




If the potential is an analytic function, we can make a Taylor expansion and write 
the integrand in the previous expression as a sum whose terms consist, in general, 
of a product between the potential and its derivatives evaluated at x and many- 
point functions, which, due to Wick's theorem, can be written in terms of two-point 
functions {y{Tk 1 )y{Tk 2 )) ■ Since these propagators are proportional to T we obtain 
an asymptotic expansion in integer powers of T whose coefficients are integrals on 
M. of products of the potential and its derivatives. 

Finally, we normalize the expectation value (1) such that for the case V(x) = 
we obtain 



(2.10) (y\e-^\ X )=e~^. { l) = e -^=- 

which is the known value for the transition amplitude for a free particle. Notice 
that the numerator in the last expression corresponds to the contribution of the 
classical trajectory whereas the denominator corresponds to the contributions of 
the quantum fluctuations. 

We end this section by recalling that for a Schrodinger differential operator 
on a manifold with boundaries the trace of the heat kernel admits the following 
small-T asymptotic expansion^ 

00 

(2.11) Tr e -n-A+v( x )} = 1 y a . T n/2 

where m is the dimension of the base manifold M.. The coefficients A n / 2 for odd n 
vanish for manifolds without boundary (see the discussion after eq. (|2.9j) .) In the 
next section it will become clear how these half-integer powers of T appear due to 
the presence of boundaries. 

3. Heat kernel expansion under Robin boundary conditions 

In this section we consider the small-T asymptotic expansion of the trace of 
the heat kernel corresponding to a Schrodinger operator on the positive half-line 



^This expansion holds under certain regularity assumptions on the manifold, the potential 
and the (local) boundary conditions. 
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A4 = R + . We impose Robin boundary conditions on the functions ip(x) in the 
domain of the Schrodinger operator 

(3.1) d x ^(x) + j-^(x)\ x=o = 
where 7 G R. 

This type of boundary conditions can be implemented in the path integral 
by computing the contributions of all paths that correspond to Neumann bound- 
ary conditions and, at the same time, introducing a Dirac delta-function in the 
potential (see appendix.) The computation of the contribution of all paths that 
correspond to Neumann boundary conditions can be performed by the method of 
images, according to which 

(3.2) (x\e- T ^ A+v ^\x) N = (x\e- T t- A+v *W\x) M + {~x\e~ T ^ A+v ^\x) R 

where the subscript N stands for "Neumann" and the subscript R indicates that 
the corresponding transition amplitudes are computed in the whole real line. Ac- 
cordingly, Vh(x) is the extension to the whole real line of the potential V(x), which 
is defined only on R + , by a reflection about the origin. As just mentioned, in order 
to obtain the transition amplitude for Robin boundary conditions we still have to 
include in the potential a Dirac delta- function. Finally, using eq. (|2.8|) to compute 
the transition amplitudes in the whole real line, the trace of the heat kernel for 
Robin boundary conditions can be written as 

Tr 7 e- T {- A+y (*» 



(3.3) 
where 

(3.4) V*(x) = V{x) + 6(-x)(V(-x) - V(x)) - 2 7 • S(x) 

The term with the Heaviside function 0(x) represents the extension of the potential 
to the negative half-line whereas the term with the Dirac delta-function is needed 
in order to obtain the transition function for the boundary conditions given by 
(|3.ip . Notice that the factor e~ x l T in the second term of eq. (|3 . 3|) is responsible 
for the appearance of half-integer powers of T. 



dx(x\e- T ^ A+v ~'^\x) R + 

OO 

dx(-x\e- T t- A+v ^ x »\x), 



°° / rp\fl POO pi pi 

V — / dx I dT X ... I dr n x 

f a n - Jo Jo Jo 



(J[VM* + y(n))j 

-e"T (f[V^ R ((l~2T k )x + y(T k )) 



\k=l 
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Next, we compute the contributions to the first coefficients A n / 2 of the heat 
kernel expansion (|2. 1 1|) due to Robin boundary conditions. Since 7 = corre- 
sponds to Neumann boundary conditions these contributions, which are propor- 
tional to positive powers of 7, must be added to the coefficients A n i 2 corresponding 
to Neumann boundary conditions. The coefficients A n / 2 for Neumann boundary 
conditions have been already computed in the worldline formalism in [Tj . 

It is interesting to notice that the inclusion of the Dirac delta-function leads 
to Dirichlet boundary conditions for 7 — > —00. Indeed, by the method of images, 
the transition function for these boundary conditions is given by the difference 
of the terms at the R.H.S. of (13. 2ft . This subtraction amounts to omit from the 
computation the contribution of all those paths that hit the origin at least once. 
In accordance, as it can be seen from the first line of eq. (|2.9p . the Dirac delta- 
function in (|3.4[) cancels the contributions of the paths that hit the origin. 

To compute the contributions of the boundary condition (|3. If) to the coef- 
ficients A n / 2 (with respect to Neumann boundary conditions) we do not use the 
procedure described in the previous section. Due to the singularities introduced in 
the potential by the Dirac delta-function and also by making the reflection with 
respect to the origin, one cannot make a Taylor expansion of V 7 ir in expression 
(|3.3[) (see the discussion after (|2 . 9[) . ) Nevertheless, one can deal with the Heavi- 
side and the delta- function introduced in (|3.4| by appropriately constraining the 
paths in the functional integration. We describe this mechanism by computing the 
contributions to the coefficients A±, A 3 / 2 and A 2 . 

The first non- vanishing contribution due to Robin boundary conditions comes 
from the term corresponding to n = 1 in eq. 
(3.5) _ 

dn {( - 2 7 • S(x + y(n)j) + e~£( - 2 7 • - 2t±)x + 2/(n)))} 

Let us consider the delta-function in the expectation value. For the first term, this 
means that one should consider the contributions of only those paths y(r) begining 
and ending at the origin in Euclidean time T such that at r = 17 reach the point 
—x. This restriction can be implemented by computing the product of the (free) 
transition amplitude (see eq. (|2.10p ) from the origin to —x in time Tti times the 
(free) transition amplitude from —x back again to the origin in time T(l — ti). The 
same calculation can be implemented for the second term in (|3.5[) but considering 
the intermediate point —(1 — 2r\)x instead of —x. The result can be written as 



r 1 

-T I dx 



2T 7 



f°° f 1 J e" 1 ^ 
/ dx dr < - 
Jo Jo I 



e 4T(1— ri) 



V^kTTl ^/At:T{1 - ri) 

(1-2t 1 )^j; z (1-2 T1 ) 2 x 2 

2 g TlWl g 4T(1- T1 ) I x 

(3.6) +e-~ e _ = / : . } = -_ ■ 2 7 



V47tTti ^47rT(l - n) V^T 
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where the result can be easily obtained interchanging the order of integration. We 
conclude that the contribution of the boundary conditions to A\ is given by 27. 

Next, we proceed with the term corresponding to n — 2 in eq. (|3.3p . There 
are two contributions to be considered, one proportional to 7 and the other one 
proportional to 7 2 . The contribution proportional to 7 2 is given by 

•4 7 2 [ dx ^ dn I dr 2 {(s(x + y(n)) ■ 5(x + y(r 2 ))) 



T 

T 



4T 2 7 2 / dx I dri / dr- 



+e-^ (S((l - 2n)x + y(n)) ■ S((l - 2t 2 )x + t/(r 2 )))} 

( x 2 

e 4TV i 1 



g 4T(1-t 2 ) 



JO Jrj 

(l-2 T1 ) 2 x 2 ( T )x 2 (1-2t 2 )^' "\ 

x 2 g «V7— e - - t" e~ 4T(l-r 2 ) T 3/2 

3.7 +e"~ ; I = -== ■ Vr?7 2 

V / 47rTri ^AttT{t 2 - n) v / 4 7 rT(l - r 2 ) J VirrT 

Therefore, the contribution of the boundary conditions to A 3 / 2 is given by y / 7T7 2 . 

The leading order in T of the contribution proportional to 7 in the n = 2 
term of eq. (|3 ,3[) contributes to the coefficient A2 and is given by 

^(-2 7 ) jf cfe^ dnjf dr 2 {(2F(x + y( n )) *(a; + y(r 2 )) 

(3.8) +e-^ (2V((1 - 2n)x + y(n)) <5((1 - 2T 2 )a; + y(r 2 )))} 

The leading order in T of this expression is obtained by making a Taylor expansion 
of the potential about the origin 

/ 2 2 

2T 2 1 V(0) 



/•<» f 1 f 1 I e -47Vi- e -4T(l-r 2 ) 

/ dx dri / dr 2 < . — 

Jo Jo Jo [ V4ttT^ ^/47rT(l - r 2 ) 



_ (1-2t 2 )^x z _ (1-2t 2 ) 2 x^ ^ 

3.9 +e~ ^ = ^ = . 27 V 

v 7 ^^??^ ^477^(1 - r 2 ) J V^T 

This is a contribution to the coefficient A 2 proportional to 7. There is another con- 
tribution to this coefficient due to the boundary conditions which is proportional 
to 7 3 . We end this section by computing this contribution, which comes from the 
n = 3 term in eq. (|3.3p and is given by 

-|^(-8 7 3 )^ dx jdn Jdr 2 Jdr 3 {(6(x + y(Ti))6{x + yfo))8(x + yfo))) 

+e-T ^((1 _ 2n )x + y(n))S((l - 2r 2 )x + y(r 2 ))S((l - 2r 3 )x + y(r 3 )))} 

T 2 4 

(3 - 10) -£rr" 
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Taking into account this last result together with the one given in eq. ()3.9|) we 
conclude that the total contribution to A2 due to the boundary conditions is given 
by 27 Vq + 4/3 7 3 . All these contributions are the expected ones (see e.g. [5].) 

4. Conclusions 

Worldline approaches can be successfully used to study QFT's on manifolds 
with boundaries. After recalling the methods employed in pQ, we have considered 
the inclusion of Robin boundary conditions. These conditions can be enforced by 
adding a delta-function potential to the setup for Neumann boundary conditions 
described in [T] . This way we have been able to recover the expected corrections to 
the heat-kernel coefficients up to, and including, A^2- We plan to further simplify 
and extend the methods reviewed in the present paper [7] . In addition it would be 
interesting to extend them to include curved boundaries, curved spaces and fields 
with nontrivial spins. 

Acknowledgments: P.A.G.P. would like to thank the participants and the or- 
ganizers of QFEXT'07. P.A.G.P. also acknowledges financial support from CON- 
ICET (PIP 6160), UNLP (11/X381), INFN and DA AD. 

Appendix A. Appendix 
Consider the following kernel, defined on the half line 
(A.l) V(Z/, x; (3) := (y\e^ H -<\x) R + (-y\e~^ \x) R 

with 

(A.2) H 7 = -d 2 y - 2 7 «%) . 

By construction, (|A.1[) is even under y — ► —y. In other words ip(y,x;/3) can be 
extended to the whole line such that it is continuous as a function of y. On the 
other hand, its first derivative is odd 

(A.3) d y ip(y, x; /3) = -d y ij}{-y, x; (5) 

and is not defined at y = because of the presence of the delta function. In fact, 
(|A.ljl satisfies the heat equation 

(A.4) (d p - 2 y - 27%)) V>(y, x;/3) = 

which is equivalent to the heat equation for a free particle away from the boundary 
whose wave function satisfies the jump condition 

(A.5) - d y ^(Q+, x; f3) + d v ^(0-,x; j3) - 2 7 ^(0, x; 0) = 

Using eqs. (|A.3j) and (|A.5|) we obtain 

(A.6) d v ^(0 + ,x;/3)+-f i/>(0,x;P) =0 

Therefore, (|A. 1|) is the heat kernel for a free particle on the half line subject 
to Robin boundary conditions. Note that taking 7^0 one correctly achieves 
Neumann boundary conditions. 
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